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ABSTRACT 

We calculate the gravitational wave signal from the growth of 10^ supermassive black holes 
(SMBH) from the remnants of Population HI stars. The assembly of these lower mass black holes 
is particularly important because observing SMBHs in this mass range is one of the primary science 
goals for the Laser Interferometer Space Antenna (LISA), a planned NASA/ESA mission to detect 
gravitational waves. We use high resolution cosmological N-body simulations to track the merger 
history of the host dark matter halos, and model the growth of the SMBHs with a semi-analytic 
approach that combines dynamical friction, gas accretion, and feedback. We find that the most 
common source in the LISA band from our volume consists of mergers between intermediate mass 
black holes and SMBHs at redshifts less than 2. This type of high mass ratio merger has not been 
widely considered in the gravitational wave community; detection and characterization of this signal 
will likely require a different technique than is used for SMBH mergers or extreme mass ratio inspirals. 
Wc find that the event rate of this new LISA source depends on prescriptions for gas accretion onto 
the black hole as well as an accurate model of the dynamics on a galaxy scale; our best estimate yields 
^ 40 sources with a signal-to-noise ratio greater than 30 occur within a volume like the Local Group 
during SMBH assembly - extrapolated over the volume of the universe yields ^ 500 observed events 
over 10 years, although the accuracy of this rate is affected by cosmic variance. 

Subject headings: galaxies, intermediate mass black holes, supermassive black holes, gravitational 
waves, dark matter halos, n-body simulations 



1. INTRODUCTION 

Any mass configuration that creates a time-dependent 
quadrupole of the stress-energy tensor will radiate 
away energy in gravitational waves. One of the 
best examples of this is the binary pulsar system 
J1915-I-1606, whose measured orbital decay of 76 /xs/yr 
is ex tremely well-described by gravitational wave emis- 
sion (|Hulse fc Tavloi|[T975f) . However, even for very mas- 
sive and compact astrophysical objects, the strain am- 
plitude of this fiuctuation in spacetime, /i, is much less 
than 10"^°. For this reason, gravitational radiation has 
yet to be directly detected by current ground-based grav- 
itational wave observatories such as the LIGO, VIRGO, 
and TAMA. 

The strongest expected gravitational wave sources oc- 
cur from the inspiral and merger of binary supermassive 
black holes. At the current epoch, nearly every galaxy is 
thought to host a supermassive black hole (SMBH) with 
a mass of 10^ - 101° M© (e.g. iKormendv fc Richstond 
|1995[ ). Observational evidence suggests that the 
black holes are embedded in galactic nuclei even at 
high redshi ft, and grow more massive as the galaxy 



ing, each galaxy m erger is e xpect e d to spawn a binary 
SMBH (iBogdanov ideFall l200l iMaver et all 120071 
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Adams et all I2001t 



■20031: IMenou et al.F 
Haehnelt fc KauffmannI 
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20041) . Since galaxies are thought to assemble by merg- 
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In general, the inspiral and coalescence of low mass 
SMBHs are expected to s pawn the large st amplitude 
gravitational wave sig nal iHaehnelti [19981 e.g.). The 
hghtest SMBH binaries, with masses lO'^ - lO"^ M©, 
have coalescence frequencies between '-^lO"* — IHz- 
squarely within the frequency range of the NASA/ESA 
Laser Interferometer Space Antenna (LISA) which has 
been proposed to launch in the next decade. LISA 
observations will complement ground-based gravita- 
tional wave observatories by broadening the observable 
gravitational wave spectrum to include low frequencies. 
In particular, gravitational wave detections made by 
LISA will be able to directly map t he assembly o f these 
lightest SMBHs (|Haehneltl [l99l IMenou et al.[ [200l 
lEnoki et al.l[200l . 

In the current picture of SMBH assembly, the black 
hole begins life as a low mass "seed" black hole at 
high redshift. It's not clear, though, when exactly 
these BH seeds emerge or what mass they have at 
birth. SMBH seeds may have been spawned from 
the accretion of low angular momentum gas in a dark 
matter halo (Kous hiappas et ahl [20041 : iBromm fc Loebl 
|2004), the coalesc ence of many s eed black holes w ithin 
a halo ( Bcgcl man fc ReesI |1978[ : [Islam et all |2004[ ). or 
from an IMBH formed, perhaps, by runaway stellar 
colhs i ons ([Portegies Zwart et al.ll2004l: IMiller fc Colbert] 
120041: Ivan der Marel[ 120041 )" However, the most likely 
candidates for SMBH seeds are the remnants that form 
from the first generation of stars sitting deep within dark 
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matter halos (iMadau fc ReesI 120011; iHeger et alj 120031: 
Volonteri et~all l2003al: llslam et al.ri2003l: iWise fc Abell 
2005f ) - so called Population III stars. With masses 
< 10^ Mq, these rehc seeds are predicted to lie near 
the centers of dark matter halos betw e en z ^ 12 — 
20 (|Bromm et all [1991 lAbel et alJl2000l [200l) . Struc- 
ture formation dictates that dark matter halos form in 
the early uniyerse and hierarchically merge into larger 
bound objects, so naturally as dark matter halos merge, 
seed black holes sink to the center through dynamical 
friction and eyentually coalesce. Dark matter halo merg- 
ers become synonymous, then, with black hole mergers 
at these masses and redshifts. This means that although 
the seed formation stops at z^l2 as Populatio n III su- 
pernovae rates drop to zero (IWise fc Abell[2005h . SMBH 
growth continues as dark matter halo mergers proceed 
to low redshifts. 

Gas accretion is thought to play a critic al role in fu- 
eling the early stages of black hole growth (Ipavid et al.l 
[19871: IKauffmann fc Haehneld[2000l : IMerlonill2004n . and 
this ma y explain the tightness of the Mbh — cr 
relati o n (iBurkert fc Silk 2001: Haehnelt & Kauffmannj 
20001: iDi Matteo et al 2005; Kazantzidis et al. 2005t 
Robertson et al. I I2006D . Since high redshift galaxies 
are thought to be especially gas-rich, each merger 
brings a fresh supply of gas to the center of the 
galaxy, and new fuel to the growing supcrmassive black 
hole (Mihos & Hernquist 1994; Di Matteo ct al. 2003). 
From a combination of gas accretion and binary black 
hole coalescence, it is thought that these Pop III- 
gene rated seeds may form the SMBHs we obserye to- 
day (ISoltanlll982l: ISchneider et al.ll2002bD . 

During a galaxy merger, each black hole sinks to 
the center of the new galaxy potential due to dynam- 
ical f riction and eventually becomes bound as a bi- 
nary (|Kai^tiidi^Fin [IqS lliciia!^^ Dy- 
namical friction then continues to shrink the orbit un- 
til the binary is hard (i.e, the separation between each 
black hole, aeBH, is such that the syste m tends to los e 
energy during stellar encounters) (Heg gie et al.l [20071 ). 
Thereafter, further decay is mediated by 3-body scat- 
tering with the ambient stellar background until the bi- 
nary becomes so close that the orbit can lose energy 
via gravitational radiation. In studies of static, spher- 
ical potentials, it may be difficult for stellar encounters 
alone to cause the binary to transition between the 3- 
body s cattering phase and the grav itational radiation 
regime (jMilosavlievic fc Merrittl[2003[ ). However, in gas- 
rich or non-spherical systems, the binary rapidly hard- 
ens and coalesces into one black hole , emitting copious 
gravitational radiation in the process (Mayer et al."2007t 
Kazantzidis ct al. 2005; Bcrczik ct al. 2006; Sigurdssoq 
20031: iHollev-Bockelmann fc SigurdssonI [20061 ). 

In our previous work, we calculated the cosmological 
merger rate for bla ck holes between 2 00 - 3 x 10^ Mq 
from redshift 49-0 ([Micic et al.l 120071 120081 (hereafter, 
M07)). Our approach combined high- resolution, small- 
volume cosmological N-body simulations with analytic 
prescriptions for the dynamics of merging black holes 
below our resolution limit; this allowed us to explore dif- 
ferent black hole growth mechanisms and seed formation 
scenarios while also accurately simulating the rich and 
varied merger history of the host dark matter halos. 



In this paper, we calculate the gravitational wave sig- 
nal from the black hole mergers involved in assembling a 
supcrmassive black hole at the center of the Milky Way 
analogue our simulation volume. The volume is designed 
to provide one possible evolutionary path for a region like 
our Local Group, and as such, it should contain super- 
massive black holes on the light end of the supcrmassive 
black hole mass spectrum - the sweet spot in SMBH 
mass for LISA observations. We include all the merg- 
ers that have occurred from redshift 49 to the present 
epoch within a 1000 Mpc'^ volume of the Universe that 
represents a Local Group type of environment. 

We found that the gravitational wave sources revealed 
in this volume are from much higher mass ratio merg- 
ers that the mergers predicted to be involved in assem- 
bling the most massive SMBHs. Most of the LISA sci- 
ence and data analysis community has been anticipat- 
ing more equal mass mergers, and have developed ex- 
tensive gravitational wave templates and parameter ex- 
traction techniques based on the assumption that the 
black hole binaries are order unity mass ratio systems 
(jBabak et al.ll20 0"8i e.g.). While this may be true for the 
most massive SMBHs, we find that the black holes in our 
volume experience mergers with mass ratios as high as 
M2/M1 = 10000 : 1. These high mass ratio mergers cer- 
tainly generate a different gravitational wave signal; they 
may even deserve a different source classification to sep- 
arate it from the classical equal mass merger or extreme 
mass ratio inspira0. 

We use the rates found in our volume to extrapo- 
late this gravitational wave signal over the Hubble Vol- 
ume. Naturally, this is highly sensitive to cosmic vari- 
ance, and ten runs selected from an initial larger volume 
are planned to mitigate cosmic variance in our merger 
rate estimates. However, these small volume cosmolog- 
ical simulations are extremely nonlinear, and as such, 
are computationally expensive. We are publishing the 
preliminary results from our first, smaller scale, simula- 
tion in order to draw the attention of the community to 
the interesting possibility that the assembly of the lowest 
mass SMBHs may involve these high mass ratio mergers 
and be a significant contributor to the LISA detectable 
event rate. 

We review the details of our simulation in section 2 
and in section 3, we describe how to calculate the grav- 
itational wave signal from two merging black holes. We 
discuss our results and implications in section 4. 

2. TRACKING THE MERGER HISTORY OF SMBHS 

In IMicic et "all ()2008l ). we performed a high- resolution 
'zoom-in' cosmological N-body simulation of a comoving 
section of a ACDM universe (f2M=0.3, f2A=0.7, as=0-9 
and h=0.7) from z— 49 to z= 0. The high resolution 
region was a box 10 h^^ Mpc on each side, for a to- 
tal comoving volume of 1000 Mpc'^. The volume is 
designed to provide one possible evolutionary path for 
a region like our Local Group, and as such, it should 
contain supcrmassive black holes on the light end of the 
supcrmassive black hole mass spectrum. 

Our mass resolution is 8.85 x 10^ M©, and our spa- 

^ Extreme mass ratio inspirals involve compact objects, like 
white dwarfs, neutron stars, or stellar mass black holes spiraling 
into a SMBH. 
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tial resolution is 2 kpc. After the simulation is com- 
plete, we identified halos with at least 32 particles us- 
ing P-Groupfinder ([Springel et al.|[200ll) and seed black 
holes those halos in the appropriate mass and redshift 
range to host Pop III stars. Note that we are using 
WMAP3 (Soergel et alJ 120071 ) cosmological parameters 
in this study to compare with our previous work; how- 
ever, at the time of this paper's submission, 'zoom-in' 
simulations of several small volumes are underway with 
WMAP5 parameters to better explore cosmic variance. 
This will allow us to better pin down the rate of these 
new gravitational wave sources. 

In our hybrid method, we combine the dark matter 
halo merger trees obtained in numerical simulations with 
an analytical treatment of the physical processes that 
arise in the dynamics of galaxy and black hole merg- 
ers. Since some of the processes are ill-constrained, we 
probe the effect of different black hole growth recipes 
on the final black hole mass function. In general, we as- 
sum e each dark matter h alo is described by an NEW pro- 
file (jNavarro et al.|[T997|) . and include the effects of dy- 
namical friction and merger-induced gas accretion onto 
the SMBH, as well as the SMBH merger itself. This hy- 
brid approach generates, for each recipe, a census of the 
number and mass ratio of the mergers in our volume at 
each redshift. 

One of the surprises from this method is that the 
black hole mass function is quite sensitive to the dynam- 
ics of the host dark matter halo on large scales. Pre- 
vious semi-analytic work on SMBH growth all assume 
that the host galaxies assembled from binary mergers 
whose timesc ale was dictated by C handrasekhar dynam- 
ical friction (jChandrasekhaii Il943f) . Numerical simula- 
tions, however, indicate that galaxies can often assemble 
from multiple near-simultaneous mergers, and that the 
merger timescale is longer than what the Chandrasekhar 
dynamical friction timescale suggests. Simulations in- 
dicate that the merger timescale is only smaller by a 
facto r 3 for minor mergers at z=0 fBovla n-Kolchin et al.l 
so given the myriad uncertainties in the dynamics 
acting in the other stages of the merger, this factor of a 
few was neglected. However, a longer dynamical friction 
timescale reduces the total number of black hole merg- 
ers in a volume of the universe by pushing some mergers 
to the future. For example, we found that the Boylan- 
Kolchin treatment for dynamical friction resulted in 1056 
SMBH mergers in our volume, while Chandrasekhar dy- 
namical friction yielded 1245 mergers. In addition, most 
of the current work on SMBH evolution indicates that 
gas accretion dominates the black hole's growth and 
that this gas is funnclcd onto the SMBH by the merger 
process itself (e.g. [R obertson et al. 2006; Hopkins et al. 
mM IDi Matteo et all 120031 : iVolonteri et al.ri2003b '): if 
the gas is driven to the center over a longer timescale, 
then each merger may excite a longer gas-fueled growth 
spurt. In the context of gravitational wave source predic- 
tion, these two effects may combine to make the SMBH 
mergers louder and at lower frequency than what has 
been predicted with a Chandrasekhar dynamical friction 
prediction. We will return to this in Section 4. 

In this paper, w e select a sample of the black hole 
growth recipes from lMicic et al.l (|2008[ ) and calculate the 
gravitational wave signal from each of the black hole 
mergers in the volume. Briefly, the recipes span two 



choices for the dynamical friction timescale, and three 
choices for the mass ratio of the halos that excite merger- 
driven gas accretion onto the black hole (4:1, 10:1, and 
oo : 1 - a.k.a. 'dry growth'). 

The next section outlines this black hole growth pre- 
scription. 

2.1. SMBH Growth Prescription 

The SMBH in our model grows through a combina- 
tion of black hole mergers and gas accretion. To bet- 
ter separate the effects of gas accretion on the black 
hole, we include a dry growth scenario, where the black 
hole grows through mergers only. At high redshift, this 
galaxy merger-driven approach is likely a good assump- 
tion, though note that at low redshift when mergers 
are infrequent, secular evolution, such as bar instabili- 
ties, may dominate the gas (and therefore black hole) 
accretion. Integrated over the whole of a black hole 
lifetime, though, this major merger-driven accretion is 
likely to be the dominant source of gas inflow. Since 
the black hole growth is so strongly dependent on what 
fuel is driven to the center during galaxy mergers, it 
is important to characterize this merger-driven gas in- 
flow, including the critical gas physics that may inhibit 
or strengthen this nuclear supply. We are motivated by 
a recent suite of numerical simulations that include ra- 
diative gas cooling, star formation, and stellar feedback 
to study the starbu rst efficiency for unequal mass ra- 
tio galaxy mergers ()Cox et al.l |2008() . which finds that 
the gas inflow depen ds strong ly on the m ass ratio of the 
galax y (see also, e.g. lHernoui st 1989; Mili os fc Hernquis^ 
[1994'). This study parametrizes the efficiency of nuclear 
star formation (i.e. gas supply and inflow), a, as a func- 
tion of galaxy mass ratio: 



0.5 



(1) 



where ao defines the mass ratio below which there is 
no enhancement of nuclear star formation (i.e. gas in- 
flow). Here, the gas accretion efficiency has a maxi- 
mum of 0.56 for 1:1 halo mergers and falls to zero at ao- 
This parametrization is insensitive to the stellar feed- 
back prescription. We use a to define how efficiently the 
merger funnels the galaxy's gas to the black hole accre- 
tio n disk. Ou r equation for a differs only slightly from 
the lCox et al.l pOOSl ) work in that we adjust ao to trigger 
gas inflow, while their ao = 0.09. 

In particular, we contrasted the black hole growth that 
would result from only major mergers with the growth 
that would occur if minor mergers were included as 
well. Our more conservative criterion, for example, al- 
lows black holes to accrete gas only if the mass ratio of 
the host dark matter halo is less than 4:1 (this is syn- 
onymous with setting ao — 0.25). This cuts off the gas 
inflow far earlier than predicted by the galaxy merger 
simulations, but is useful as a conservative estimate and 
to compare to our previous work. 

We can set an upper constraint on the final black hole 
mass by allowing gas accretion as long as the merging 
dark matter halos have a mass ratio less than 10:1 (con- 
sistent with ao = 0.1). Note that ao = 0.09 would im- 
ply that mass ratios less than roughly 11:1 would drive 
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gas inflow ()Cox et al.ll2008D . rather than 10:1. Adopt- 
ing ao = 0.09 would increase the black hole mass only 
slightly, and since our main goal is to compare to our 
earlier work, we retain ao = 0.1. 

We used a semi-analytic formalism to calculate the 
dynamical friction decay time and subsequent merger 
timescale of each black hole once it enters a dark mat- 
ter halo. Now that we have a realistic description of the 
merger time for each black hole within a halo, we can 
allow the black holes to grow for a physically-motivated 
accretion timescale. The accretion of gas onto both in- 
coming and the central black hole starts when the two 
black holes are still widely separated, at the moment of 
the first peric enter passage, and continues until the black 
holes merge (jPi Matteo et all l2005t iColpi et al.l 120071 
c.f.). This sets the accretion timescale, tacc, as follows: 
tacc = tdf(r=Rvir) - tdyn(r=Rvir), whcrc tdf(r=Rvir) is 
the merger timescale including dynamical friction; and 
tdyn(r=R.vir) IS dynamical time at virial radius Rvir, 
which marks the first pericenter pass of the black hole. 
By stopping the accretion as the black holes merge, we 
roughly model the effect of black hole feedback in stop- 
ping further accretion. 

Putting these pieces together, the mass accreted by a 
black hole during tacc(i'=Rvir) is: 
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where Mb h,o is initial black hole mass, and a is starburst 
efficiency (jCox et al.l [20081 ) . and tsai = eMBH,ocV[(l - 
e)L], where e is the radiative efficiency, L is the luminos- 
ity, and c is the speed of light. After tdf , the incoming 
black hole merges with the SMBH at the center and a 
new SMBH is formed after having accreted gas for tacc- 
The accretion time and efficiency both implicitly encode 
the large-scale dynamics of the merger and the bulk gas 
accretion into the nuclear region, while tgai describes the 
accretion disk physics. As before, we set tsai to describe 
sustained Eddington-limited accretio n with a efficiency 
of e = 0.1 (IShakura fc Sunvaevl [l973h. 

In lMicic et al.l (I2008D . we found that a gas accretion 
triggered by major mergers (4:1) and a Boylan-Kolchin 
dynamical friction prescription produced a 4 x 10^ Mq 
SMBH from a 200 Mq seed in place by z= 5. We con- 
sider this our preferred model as it produces a SMBH 
consistent with the Sgr A*. Note, however, that the mi- 
nor merger (10:1) Boylan-Kolchin prescription also gen- 
erates a 1.3 X 10^ M0 black hole from the same seed, and 
can be considered a viable model for an M31-like SMBH. 

3. CALCULATING THE GRAVITATIONAL WAVE SIGNAL 

In order to compare our results with other published 
results, we follow the approach of ISesana et al, (20041 
to calculate the gravitational wave signal from the black 
hole mergers in our simulation volume. In this section, 
we provide a brief background on the theory and outline 
the method; for a more in-depth description of gravita- 
tional waves from SMBHs see, for example, ^Berti (2006). 

The orbital motion of the binary system constitutes 
will excite a time-dependent mass quadrupole, which 
generates gravitational radiation. The gravitational 
wave strain, h, measures the strength of the propagating 
wave. The amplitude depends on the comoving distance 



to the binary, d, the rest frame gravitational wave fre- 
quency, fr, and the chirp mass of the binary, A4. The 
rest-frame gravitational wave frequency is related to the 
orbital period, P, of the binary, fr = 2/P, and the chirp 
mass depends on the mass of each binary component as 

follows: Ai = (TOim2)'^/'^/(mi -I- 7712)^/^. When averaged 
over the sky position, polarization, and period, the strain 
can be written as: 
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(3) 



This gravitational radiation will cause the binary com- 
ponents to inspiral on ever tighter and faster orbits. 
When massive black hole binaries are widely separated 
and orbiting at low frequencies, the gravitational radi- 
ation emitted is relatively weak, and therefore the fre- 
quency shift per orbit is minuscule. Most of the evolution 
is spent in this phase, where a frequency shift of order 
unity takes many orbits to achieve. 

As the binary separation slowly shrinks, the gravita- 
tional radiation emitted increases dramatically until the 
binary coalesces. Close to coalescence, the binary rapidly 
sweeps through many frequencies in one orbit, and at the 
innermost stable circular orbit (ISCO), the change in fre- 
quency per orbit is of order unity. 

We set the minimum observed frequency to /isco for a 
test particle orbiting around a single Schwarzschild black 
hole. In the black hole mergers we consider, the mass 
ratios are always well outside this test particle limit; we 
adopt the conventional definition for /isco nonetheless: 



/isco 
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where c and G are the speed of light and gravitational 
constant, respectively, mi and 1112 are the black hole 
masses, and z is the redshift. 

The frequency shift rate in the rest-frame is, to first 
order: 

fr = /,"/3_A/(5/3_ 

The orbital periods of million solar mass black hole bina- 
ries in the gravitational radiation regime are on the order 
of hours, which means that multi-year observations could 
accumulate many cycles of an inspiral at a particular fre- 
quency. LISA, for example, is expected to observe the 
sky for at least a year, and projections place the LISA 
lifespan at roughly a decade. Million solar mass black 
holes orbiting in the LISA frequency band of ^ 10"'' — 1 
Hz are years or less from coalescence, so LISA observa- 
tions will be able to track the inspiral and coalescence 
phases of supermassive black hole binaries. The number 
of cycles accumulated during the inspiral in an observa- 
tion of duration r depends on the observed gravitational 
wave frequency, n = fr, where / = /r/(l + z). Note 
that there are two regimes: / > n/r, where the binary 
sweeps through many frequencies in an observation, and 
/ < n/r, where the signal builds from several orbits at 
one frequency. 

The observed characteristic strain, he, is the strain ac- 
cumulated a single observation: 
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Since we have the component mass and redshift of each 
black hole merger in our simulation volume, we can cal- 
culate h and he over an observation span r before merger. 
Most black holes in our volume will merge within a year 
of reaching the LISA band, so a 3-year observation win- 
dow should catch these mergers in the act - if this volume 
is a representative slice of the Universe. We distinguish 
between the mergers directly found in our volume and 
the mergers extrapolated throughout the Universe in Ta- 
ble H 

4. RESULTS 

4.1. A New Gravitational Wave Source 

Overall, most of the mergers in our volume are between 
seed and intermediate mass black holes (O(IO^)Mq) at 
redshifts greater than 5 (Figure [T]). Though these merg- 
ers were critical in assembling the eventual SMBH in our 
Milky Way analogue, the generated LISA signal-to-noise 
ratio of most of these low mass mergers at high redshift 
were much less than 1.0. Figure [7] demonstrates that 
the overwhelming majority of black hole mergers will fall 
below LISA'S detection limit. Although there are a few 
mergers that are resolvable by LISA at z> 6, the number 
of resolvable sources increases once the SMBH grows to 
a mass within the LISA band at z < 5, as can be seen i n 
Figure [5] This broadly agrees with lSesana et all (|2005[ ). 
which showed that ~ 10% of the precursors to 10^ M© 
SMBHs are expected to be observable with LISA out to 
redshift 10. 

However, Figure [T] also reveals a substantial class of 
black hole mergers with very high mass ratios; these arise 
from the accretion of smaller satellite halos at low red- 
shift. Since these mergers are nearby, they are easily 
detectable with LISA. This is in stark contrast to the 
predictions for assembling the mass ive end of the SMBH 
mass spectrum; iSesana et al.l (|2007f) showed that the de- 
tectable mass ratios are equally distributed in the range 
of 1-10. In fact. Figure m indicates that the most com- 
monly observed black hole mergers in our volume have 
mass ratios of 1000 or more (depending on the black hole 
growth prescription). Mergers of more equal mass dark 
matter halos (and subsequently, the coalescence of more 
equal mass black holes) occurred in this volume at 2; > 8, 
and at this epoch, most of our black holes were too small 
to be observed by LISA. Figures [T] and [3] bear this out 
by plotting mapping the number of resolvable sources as 
a function of mass ratio and redshift for one black hole 
growth prescription over a 10 year observation. 

The signals from an equal mass inspiraling system and 
those with mass ratios of upwards of ten thousand can 
be quite different. From an analytical standpoint, the 
difference in signals can be understood by investigating 
the post-Newtonian expansion of the gravitational wave 
strain, which uses the orbital velocity as an expansion 
parameter. Each term in the expansion results in har- 
monics of the orbital frequency .Blanchet et al.l (pL996). 
Starting at the first full order, certain select frequen- 
cies are scaled by 77 = mim2/(mi + m-i)^ , which sup- 
presses those frequencies for systems with large mass ra- 
tios. Figure [6] compares the power spectral densities for 
two SMBH systems: one with two lO"^ Mq SMBHs and 



the other with mi = 10^ Mq and m-i = 10'* Mq. These 
figures demonstrate how the equal mass system sweeps 
through all frequencies while the large mass ratio system 
shows much more structure, with many of the frequencies 
being suppressed. 

The LISA data streams present unique challenges for 
data mining. Unlike electromagnetic observatories, LISA 
is simultaneously senstitive to gravitational wave sources 
located all throughout the sky. LISA will also be sensi- 
tive to a wide variety of astrophysical sources - from su- 
permassive black hole binaries, to millions of close white 
dwarf binaries within our galaxy, to EMRIs. The data 
analysis objective is to conclusively detect a signal and 
to extract descriptive parameter values. In preparation 
for the immense data analysis challenge, simulated data 
has been produced and distributed to the LISA commu- 
nity as pa rt of the Mock LISA Data Challenge (MLDC) 
iVallisneril (i200_9)- So far the challenges have focused on 
order unity mass ratios for the simulated supermassive 
black hole inspirals. While for most analysis techniques 
an assumption about the mass ratio is not hardwired into 
the routine, the routines have not been tested at the more 
extreme ratios of 10^ as suggested here. It is possible 
that the surpressed harmonics will cause false-positives 
or negatives in some routines because the system will be 
missed or misidentified. 

4.2. The LISA Signal from Assembling a Milky Way 

SMBH 

Figure |9] presents the total characteristic strain for 7 
bright black hole mergers in our volume, assuming a 10:1 
black hole growth recipe and a Boylan-Kolchin dynami- 
cal friction treatment. The differences in the two classes 
of source is clear - the brightest sources are high mass 
ratio mergers which coalesce in the LISA band at low 
redshift, while the other bright class probes more equal- 
mass mergers of low mass black holes at high redshift. 
Note that the black holes in the bottom-most track actu- 
ally merge outside the LISA band and will be considered 
an inspiral; this was the only detectable inspiral in our 
volume. 

In Table [U we determined the redshift at which these 
7 mergers would become undetectable in a 3 year obser- 
vation. Understandably, the highly unequal mass merg- 
ers can only be detected out to a luminosity distance of 
roughly 3 Gpc, which is similar to the distance probed 
by extreme mass ratio inspirals (|Gair "erall[200l . LISA 
observations of this class of high mass ratio merger, then, 
may be useful to probe the fraction of black hole mass is 
accreted at late times for this low mass SMBH range. 

Figure [5] presents the total characteristic strain for the 
LISA detector from all the black hole mergers in our 
1000 h~3 Mpc'^ volume over a 3-year observation span. 
Each curve has a low frequency rise that drops steeply 
off before hitting a shallow plateau; the rise and drop- 
off is the signature of the accumulated high mass ratio 
mergers, while the shallow segment marks the sum of 
the relatively equal mass mergers. We caution that this 
implicitly assumes that all of these mergers would take 
place during this 3-year period, or alternatively that this 
volume is representative of the Universe. Issues of cosmic 
variance aside, we can see that the black hole growth 
prescription strongly influences the total strain in the 
LISA band; this is best seen when dynamical friction 
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is described by iBovlan-Kolchin et alj ()2008l ) where the 
drop-off changes by a decade in frequency in response to 
a change in the typical mass of the local SMBHs. 

4.3. Extrapolating to Universal Black Hole Merger 

Rates 

Although this is plagued by cosmic variance, it is in- 
structive to estimate the number of mergers observable 
by LISA in the Universe expected from assembling these 
lightest SMBHs. Table [5] presents the extrapolated Uni- 
versal LISA merger rates. One surprising point is that 
the LISA merger rates depend so strongly on the adopted 
form of the dynamical friction force. All previous LISA 
merger rate estimates have used a semi-analytic tech- 
nique that employs Chandrasekhar dynamical friction 
to merge the dark matter halos and to usher the black 
holes to the i nner kiloparsec. However, both pertur- 
bation theory (iColpi et a l.lll999 D and numerical simula- 
tions (Wcinbcra" 19891 iH oUcv- Bockelmann fc Richstond 
119991 [Bovlan-Kolchin et al.ii2008i ) have shown that Chan- 
drasekhar dynamical friction approximates the merger 
time to within a factor of two, at best. When we em- 
ploy a dynamical friction formalism that is based on fits 
to merger timescales from numerical simulations, wc find 
that the black hole mergers are delayed to lower redshift. 
Naively, this would simply make every merger louder. 
However, in our SMBH growth prescription, there is an 
additional effect: since the incoming SMBH drives gas 
infiow to the primary SMBH over a longer timespan, re- 
sulting SMBH ultimately grows more massive than with 
Chandrasekhar dynamical friction. This can place the 
SMBHs in our volume just out of LISA's 'sweet spot' 
which will decrease the signal-to-noise ratio of the more 
massive black holes. Paradoxically, then, the LISA event 
rate for more correct merger times drops by as much as 
an order of magnitude. For our most realistic model to 
assemble the lightest SMBHs - gas accretion triggered by 
major mergers, with a Boylan-Kolchin dynamical friction 
merger timescale - the merger rate drops by a factor of 
two. Over a 3 year LISA observation, we should be able 
to detect at least 70 SMBH mergers with a signal-to-noise 
ratio greater than 30 that are involved in assembling the 
light end of the SMBH mass spectrum in the Universe. If 
we consider longer LISA observation windows, the num- 
ber of observed sources increases, naturally, and for a 10 
year observation, LISA should detect nearly 500 mergers 
from the assembly of the lightest SMBHs alone. 

5. DISCUSSION 

By concentrating on a small cosmological volume, we 
have been able to model the gravitational wave sources 
that result from the assembly of a '-^ 10^ — 10^ M0 black 
hole in a Milky Way mass halo. We have calculated the 
gravitational wave strain for each of the black hole merg- 
ers in our volume, and have determined which of those 
will be detectable with LISA. Of the 1500 mergers in our 
volume, we uncovered approximately 300 — 1200 merg- 
ers detectable with a signal-to-noise ratio greater than 5 
over a LISA observation span of 3 years. 

We found that the most common class of observable 
black hole merger in our volume is between a SMBH and 
IMBH (of 200-2000 Mq ) at z< 0.05. These IMBHs origi- 
nally resided in small dark matter halos that merged with 
the massive primary halo at high redshift and had very 



long dynamical friction timescales. Before the merger oc- 
curs, the incoming IMBH may be observed with the next 
generation of X-ray telescopes as a ULX source with a 
rate of about ^3-7 yr~^ for 1 < z < 5. Because 
of their potential tie to observable ULXs, and because 
this class of source has a different waveform character 
than other well known gravitational wave sources, such 
as equal mass mergers, intermediate mass ratio inspirals 
(IMRIs), or extreme mass ratio inspirals (EMRIs), we 
have nominally dubbed this class of source as an Ultra 
Large Inspirals (ULIs). The other class was IMBH- IMBH 
mergers at z= 2 — 8. 

Note that we could not resolve the growth of the most 
massive dark matter halos with our technique, and miss 
the mergers that arise from the assembly of the most 
massive SMBHs. Therefore, we co nsider our approach 
to be complementary to studies like ISesana et al.l (j2005l 
[2OOI - we believe that the high mass ratio mergers iden- 
tified here would add to the rates found in previous stud- 
ies that focus on black hole growth in present-day halos 
larger than ^ 10^^ Mq. However, the high mass-ratio 
channel that we have identified may well dominate the 
LISA black hole merger events. In fact, even our most 
pessimistic estimate (~ 70 events) yields a com parable 
number of LISA sources as ISesana et al.l ((20051) (~ 90 
events) . 

The signal-to- noise ratio listed in Table [1] assumes sim- 
ple data analysis techniques. If we were to employ similar 
sophisticated tools as are being developed for the extreme 
mass ratio inspirals, such as matched filtering, it is pos- 
sible that this class of source may probe slightly larger 
distances than EMRIs. In any event, current data anal- 
ysis techniques to extract SMBH signals from LISA data 
streams all assume the binaries have mass ratios close to 
unity. As we have shown here, the more complex wave- 
form structure of these ULIs may require a different data 
analysis strategy. 

In our volume, the SMBH was in place at nearly its 
final mass at redshift 5, and it was built by mergers of 
hundreds of 0(10^ ) Mq black holes at z> 5 - these merg- 
ers were too weak to be detectable by LISA. In fact, the 
first trace of the growing SMBH occurs at about redshift 
7 for our most aggressive gas accretion prescription. This 
may have implications for how well LISA observations 
can constrain the early growth of these lightest SMBHs. 

In this paper, we have neglected gravitational wave 
recoil, a potentially important mechanism that may in- 
hibit black hole growth. Binary black holes strongly 
radiate linear momentum in the form of gravitational 
waves during the plunge phase of the inspiral - result- 
ing in a "kick" to the new black hole. This, in itself, 
has long been predicted as a consequence of an asym- 
metry in the binary orbit or spin configuration. Previ- 
ous kick velocity estimates, though, were either highly 
uncertain or suggested that the resulting gravitational 
wave recoil velocity was relatively small, astrophysically 
speaking. Now, recent results indicate the recoil can 
drive a gravitatio nal wave kic k velocity as fast as ^ 4000 
kms~^ (e.g. Herrman n et a l. 2007; Gonzalez ct al. 200<i 
Gonzalez et al. 2007; Baker et al. 2006; Koppitz ct aQ 
2007c iCampanelli et al.i i2007i : iSchnittman fc Buonannol 
2007f ). In reality, much smaller values than this max- 



imum may be expected in gas-rich galaxies due to the 
alignment of the orbital angular momentum and spins 
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of both black holes (|Bogdanovic et al]|2007f) . However, 
even typical kick velocities (~ 200 kms""'^) are inter- 
estingly large when compared to the escape velocity of 
typical astronomical systems - low mass galaxies, as an 
example, have an escape veloc ity of ~ 200 kms~^ (e.g. 
iHollev-Bockelmann et al.]|2008D . The effect of large kicks 
combined with low escape velocity from the centers of 
small dark matter halos at high redshift plays a major 
role in suppressing the growth of black hole seeds into 
SMBH. Even the most massive dark matter halo at z>ll 
can not retain a black hole that receives > 150 kms~^ 



kick (iMerritt et all [200l iMicic et al.l 120061 ). We have 
submitted a companion paper that incorporates the ef- 
fect of recoil velocity on the expected merger rates and 
the growth of a Milky Way-mass SMBH, and find that 
if there is no spin alignment mechanism, then a Pop HI 
seed black hole can reach 10^ M© only 20% of the time 
through merger-driven gas accretion. We are exploring 
the effect of recoil on the gravitational wave signal in a 
forthcoming paper. 
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Fig. 1. — Map of the number of black hole mergers in our 1000 Mpc^ volume as a function of the log of the black hole mass ratio and 
the redshift of coalescence. The colors represent the number of black hole mergers, with a scaling as denoted by the colorbar to the left of 
the map. Here, the dynamical friction is treated with a Boylan-Kolchin formula and gas accretion is triggered for dark matter halo msss 
mergers with mass ratios smaller than 4:1 - major mergers. Note the bimodal distribution, and the fact that most mergers are low mass 
and large mass ratio. 
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Fig. 2. — Map of the log of the signal-to- noise ratio as a function of the log of the black hole mass ratio and the redshift of coalescence. 
The colorbar to the left translates the colorscale to the log of the SNR - an SNR of 5 would be dark blue. Here, the dynamical friction 
is treated with a Boylan-Kolchin formula and gas accretion is triggered for dark matter halo mergers with mass ratios smaller than 4:1 - 
major mergers. Over this 10 year LISA observation, we can expect to detect no equal mass mergers, but that the loud sources will all have 
mass ratios larger than 100. We reemphasize that this map tracks the sources only from the formation channel of the lightest SMBHs. 
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Fig. 3. — Map of the number of resolvable sources as a function of the log of the black hole mass ratio and the redshift of coalescence. 
The colorbar to the left translates the color to the number of resolvable sources with a signal-to-noise ratio is greater than 5. Here, the 
dynamical friction is treated with a Boylan-Kolchin formula and gas accretion is triggered for dark matter halo mergers with mass ratios 
smaller than 10:1 - minor mergers. Over a 10 year LISA observation, this model is dominated by high mass ratio mergers. 
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Fig. 4. — Histogram of the mass ratios of sources predicted with a signal-to-noise ratio larger than 5, assuming a 3 year LISA observation. 
The top panel assumes that the dynamical friction timescale is set by the Chandrasekhar dynamical friction formula, while the bottom 
panel assumes the Boylan-Kolchin formula. The colors correspond to different black hole growth assumptions: red assumes black holes 
grow only through mergers, black assumes that gas accretion onto both black holes is triggered if the host dark matter halos have a mass 
ratio smaller than 10:1, and the blue assumes this accretion only occurs for halo mergers with 4:1 mass ratios. For more detail on the 
particular black hole growth prescription, see Micic, HoUey-Bockelmann, and Sigurdsson 2008. Regardless of the growth prescription or 
the large-scale dynamics, we see here that LISA observations will be dominated in number by high mass ratio black hole mergers, at least 
for the assembly of these relatively light supermassive black holes. 
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Fig. 5. — Total characteristic strain as a function of observed frequency for all the massive black hole mergers in the cosmological volume 
for a three year observation. The solid curve is the sensitivity curve for LISA for a single-arm Michaelson configuration. 
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Fig. 6. — The plus polarization of the gravitational wave amplitude as a function of frequency for two representative black hole mergers 
at redshift 3. Sources are placed at the same random sky position and orientation. Left: merger with black hole masses Ml = 1O^M0 and 
M2 = IO^Mq. Right: merger with equal black hole masses of 10^A/q. Notice the broader bandwidth and richer structure of the high mass 
ratio merger. 
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Fig. 7. — Histogram of the signal-to-noise ratio of the black hole mergers in our volume (scaled to cosmological volumes) for an assumed 
3-year LISA observation window. Sec figure |4] for a description of the panels and colors. Since our cosmological volume is small, we are 
unable to resolve the growth of the most massive black holes. Most of our mergers are at high redshift between seed and intermediate mass 
black holes involved in assembling the light end of the supermassive black hole spectrum; these high redshift, low mass mergers yield low 
signal-to-noise sources. 



TABLE 1 

Maximum redshift to observe mergers at a signal to noise ratio of 5 in a three year LISA observation. 



^^bhi ^^bh2 redshift in volume max redshift 



19398 


16654 


8.881 


9.1 


83201 


20548 


7.958 


17.5 


2723 


929 


0.726 


0.8 


46124 


9093 


1.169 


9.8 


13054368 


435 


0.017 


0.5 


13055368 


200 


0.006 


0.4 


13055168 


600 


0.006 


0.5 



Column 1: Primary black hole mass. Column 2: Secondary black hole mass. 
Column 3: Redshift the merger occurred within the simulation volume. Col- 
umn 4: The redshift to which the merger can be detected at a signal to noise 
ratio of 5 for a three year LISA observation. 
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Fig. 8. — Histogram of the redshift of the resolved black hole mergers in our volume (scaled to cosmological volumes) for an assumed 
3- year LISA observation window. See figure |4] for a description of the panels and colors. We assume that a signal-to-noise ratio greater 
than 5 will be resolvable; given the higher mass ratios of a typical merger, however, signal-to-noise ratios higher than 30 may be needed to 
resolve a source. 
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Fig. 9. — Characteristic strain over a 3- year observation for selected classes of resolvable black hole mergers in the simulation as a function 
of observed frequency. Here, the dynamical friction is treated with a Boylan-Kolchin formula and gas accretion is triggered for dark matter 
halo mergers with mass ratios smaller than 10:1. The merger redshift and pre-merger masses of the black hole binary are labeled; our 
loudest sources are either caused by mergers between two intermediate mass black holes, or between a very light seed black hole and the 
most massive supermassive black hole in our volume. 
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TABLE 2 

Number of Resolvable Black Hole Mergers For Several Black Hole Growth Scenarios. 
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Columns 1 and 2 describe the blax;k hole growth model; column 1 indicates which halo mass ratio would trigger gas 
accretion onto the primary bla<;k hole; column 2 indicates the dynamical friction formalism that dictates the gas 
accretion timecscale. Column 3: Signal to noise ratio for detection. Column 4: Observation duration. Column 5: 
Predicted number of LISA events. Note that column 5 is not the entire LISA event rate; it is the conjectured rate from 
the class of sources involved in assembling the lightest SMBHs. Also note that although the predicted event rate varies 
here by a factor of 50, even the most pessimistic rate estimate still will provide ample new sources for LISA to observe. 



